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Abstract 



Using the values of the cosmological parameters obtained from Wilkin- 
son Microwave Anisotropy Probe (WMAP) and the cosmological chemical 
clock observed by high-redshift quasars, it is shown that star formation 
should have started by the redshift zp > 9.5, and that zp is sensitive only 
to rimo^^ in flat ACDM cosmology. 
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The first year data release of Wilkinson Microwave Anisotropy Probe (WMAP) has 
opened a new era of high precision cosmology (Bennet et al. 2003). As for the values of 
the cosmological parameters it has basically confirmed that of the concordance model with 
significantly smaller error bars than before. Together with the observation of nearly scale- 
invariant adiabatic fluctuations, the inflation-based flat ACDM model can be regarded as 
the primary candidate of the standard evolution model of our Universe. 

Ever since Bubble's discovery of cosmic expansion (Hubble 1929), who concluded 
that the cosmic expansion rate was as large as Hq ~ 500km/s/Mpc, we encountered 
the cosmic age problem from time to time in the history of modern cosmology. The 
introduction of vacuum-like energy density with negative pressure such as cosmologi- 
cal constant A as inferred from the magnitude-redshift relation of high-redshift type la 
supernovae (Perlmutter et al. 1999; Ricss et al. 1998) has reconciled a relatively large 
Hubble parameter h = iJo/(100km/s/Mpc)^ 0.7 (Freedmann ct al. 2001) with a fairly 
large cosmic age to ^ 13Gyr. Indeed the cosmic age obtained by WMAP observations 
supplemented by other cosmological observations, to — 13.7 ± 0.2Gyr with h — 0.7llo]o3) 
^mo^^ = 0.135^0 009, fivo = ^toto ~ ^^mO; ^nd iltoto — 1-02 ± 0.02 are in good agreement 
with astrophysical estimate of globular cluster age and cosmological nuclear chronology 
(Spergel ct al. 2003). Here Q^o denotes the current value of the vacuum-like energy 
density or dark energy in unit of the critical density. 

Thus as far as the current Universe is concerned we have indeed the concordance of 
the cosmic age and various cosmological parameters at hand. However, the situation may 
change drasitically if we turn our attention to the high-redshift Universe. Indeed with 
the improved accuracy of current values of the cosmological parameters, we can calculate 
the cosmic age at redshift z, t{z), more accurately. Comparing t{z) with other measures 
of cosmic age at high-redshift we can obtain important cosmological and astrophysical 
information and constraints. 

As for the latter we consider observation of Fc/Mg abundance ratio probed by high- 
redshift quasars. This quantity can serve as a cosmic chemical clock due to the difference 
of the origin (Hamann & Ferland 1993). That is, Fe is predominantly produced by type 
la supernovae whose precursor has a lifetime of ~ 1 Gyr or larger, while Mg and other 
alpha elements are produced from type H supernovae with much shorter time scale. Thus 
the ratio Fe/Mg at the redshift z indicates the time elapsed between the star formation 
epoch zp and z. 

Previously this measure was mainly used to constrain the values of cosmological param- 
eters with particular emphasis on whether cosmological constant is nonvanishing (Yoshii, 
Tsujimoto, & Kawara 1998). But at the new era of precision cosmology with WMAP 
and other sophisticated observational tools we can extract useful information on zp and 
chemical evolution models of the Universe gued below. 

Although the observational constraint of WMAP on the spatial curvature or the total 
energy density of the Universe still has a two percent uncertainty, fltoto — 1-02 ± 0.02, 
we assume Jltoto is equal to unity with much higher accuracy, because we have good 
reasons to believe in standard inflation as argued above, which predicts Ifitoto ~ 1| < 
Q — 5 X 10~^ with Q being the quadrupole anisotropy of cosmic microwave background 
radiation (Kashlinsky, Tkachev & Frieman 1994). Then, in the spatially flat Robsertson- 
Walker spacetime with the scale factor a{t), cosmic age at redshift z is given by 
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r da 1 r°° dz 
^ J ~h ^ irjz {l + z)E{z)' 

r 1 1/2 

E(z)^[n^o(i + z)' + nyof{z)] , (1) 

where f{z) is a function which specifies the evolution of dark energy component. If it has 
a constant index of equation of state, Wy = Pv/Pv, where p and p denote energy density 
and pressure respectively, one finds f{z) — {1 + z)^^^^'^'"\ The pure cosmological constant 
corresponds to -u^t, = — 1. Below we mainly consider this case because it is theoretically 
well-motivated (Yokoyama 2002a,b) and observationally supported among many other 
candidates of dark energy (Yokoyama 2003). In the above expression we have neglected 
radiation component, which is a good approximation in the range of redshifts we consider. 

Several authors have calculated time evolution of Fe/Mg ratio using a chemical evo- 
lution model. For example, Yoshii, Tsujimoto & Kawara (1998) calculated its evolution 
based on the model of Yoshii, Tsujimoto & Nomoto (1996) with the following three basic 
ingredients. The first one is the star formation rate C{t) parametrized as C(t) = [fg{t)]^ 
where fg(t) is the gas fraction and and k are parameters. As the standard model 
they choose k = 1 and z/fc=i = 7.6Gyr~^ as inferred from metal abundance in quasar 
host galaxies (Hamann & Ferland 1993). The second is the initial stellar mass function 
(IMF) which is assumed to be time-invariant with a spectrum (f){m)dm oc m~^dm for 
O.OSMq < m < 50Mq with x = 1.35 based on Tsujimoto et al. (1997). The third and 
the most important ingredient is the fraction. A, and the lifetime, tja, of the progeni- 
tors that eventually produce SNe la, and their time spread. They adopt A = 0.055 for 
m = 3 — SMq and A = outside this mass range as a standard value, tja is obtained from 
the break in [0/Fc] observed at [Fe/II]~ —1 in the solar neighborhood (Yoshii, Tsujimoto 
& Nomoto 1996) and taken as tia = 1.5 Gyr, while its spread function, g{tia), is modeled 
as a power-law g(tia) oc t]^. In practice, they take 7 = and g{tia) 7^ for = 1 — 3 
Gyr, again to match the observational features of [0/Fe] and [Fe/H]. 

As for the observation, Yoshii, Tsujimoto and Kawara used UV-optical spectrum of 
QSO B1422+231 at z = 3.62 and found FeII(UV+opt)/MgII AA 2796, 2804 flux ratio 
to be 12.2 ± 3.9. Assuming this ratio traces the abundance ratio Fe/Mg, they obtained 
[Mg/Fc]= -0.61 

-a30-ai2) which, according to their standard evolution model, implies 
that the time elapsed from zp to z — 3.62 is given by 

At^=3.62 = t{z = 3.62) - t{zF) = 1.50 Gyr. (2) 

They have also considered possible modification of their standard model and concluded 
that the lower bound is 1.3 Gyr, namely, Atz=3.62 > 1-3 Gyr, which may be realized for a 
higher value oi A> 0.08. As is seen in Table 1 of Yoshn, Tsujimoto and Kawara (1998), 
modification of other model parameters tends to increase At. 

This quasar is by no means the only QSO whose flux ratio FcII/Mgll has been ob- 
served. Thompson, Hall and Elston (1999) analyzed spectra of 12 high-redshift quasars 
to probe evolution of Fell/Mgll ratio. As a result they found no decline in the iron 
abundance even at 2; = 4.47 and concluded that about IGyr had passed from zp to 
z — 4.47. Dietrich et al. (2002), on the other hand, observed six high-redshift QSOs with 
z = 3310 — 3.488 and basically confirmed the conclusion of Yoshii, Tsujimoto and Kawara 
(1998), namely, At^^^sA ~ 1-5 Gyr. 
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Comparing these values with the formula (1) calculated with the improved values of 
cosmological parameters obtained with WMAP, we can find constraint on the epoch of 
star formation zp- Specifically we consider the three possible tests on zp as argued above, 
to draw contours of zp corresponding to /S.tz=zs2 = 1-3 Gyr, Atz=3S2 = 1-5 Gyr and 
Ai,=4.77 = 1-0 Gyr. 

Figures 1-3 depict these contour on an Q^q — Hq plane. As is seen there the condition 
Atz=3S2 = 1-5 Gyr is more stringent than that Atz=4,77 = 1.0 Gyr. For the most likely 
value of the cosmological parameters obtained by WMAP, h = 0.71 and Oj„o = 0.268 
under the assumption that fl^o = 1 — ^rao, we find zp = 10.13 from Atz=3.62 = 1-3 Gyr, 
zp = 14.98 from Atz=3.62 = 1-5 Gyr, and zp = 11.91 from Ai^=4.77 = 1-0 Gyr. Thus 
significant star formation before zp > 10 is inferred. In principle we could yield even 
more stringent constraint than zp = 14.98, if we applied Yoshii et al's model (1998) to 
Thompson et al's data (1999), but we remain conservative here, allowing the room for 
uncertainty in chemical evolution models. 

More interesting are contour maps depicted on a ilmo^^ ~ Hq plane, which are shown 
in figs. 4-6. As is seen there in the range of cosmological parameters under consideration 
Zp is insensitive to Hq and depends only on Q^^^Qh"^. Thus for fimo^^ — 0.135tooo9, we 
obtain zp = 10.lio;g from Atz=3.62 = 1-3 Gyr, zp = 15. Ol];;^ from Atz=3.62 = 1-5 Gyr, and 
Zp = 11.9 ± 0.7 from At^=4.77 = 1-0 Gyr. We may therefore conclude zp > 9.5 in the fiat 
ACDM Universe with the cosmological parameters obtained by combination of WMAP 
and other measures. 

Finally we briefly mention the case dark energy has a softer equation of state w > —1 
than the cosmological constant. For a constant coefficient w, WMAP has obtained a 
constraint w < —0.78 (Spergel et al. 2003). Dark energy with a softer equation of state 
tends to require star formation at higher redshift, although the change is modest. For 
example, for w — —0.8, h = 0.71, and Q^o — 0.268 we find zp = 10.36, 15.57, and 
12.10, respectively instead of zp — 10.13, 14.98, and 11.91 for the case of the cosmological 
constant as discussed above. 

In conclusion we have pointed out that in the new era of precision cosmology we 
can obtain some useful information on the evolution of the high-redshift Universe from 
the cosmological chemical clock rather than constraining current values of cosmological 
parameters with them. We have shown that the epoch of star formation zp should satisfy 
Zp > 9.5 and that it is sensitive only to Q^ioh'^ in the context of the standard flat ACDM 
cosmology. This constraint would bring about interesting implication to formation of 
galaxies that host quasars. 

This work was partially supported by the JSPS Grant-in-Aid for Scientific Research 
No. 13640285. 
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FIG. 1. Contour map of zp which satisfies Atz=3.62 = 1-3 Gyr as a function of fl^Q and Hq. 




FIG. 2. Contour map of zp which satisfies /S.tz=3.62 = 1.5 Gyr as a function of OmO and Hq. 
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Figure 3 

FIG. 3. Contour map of zp which satisfies Ai2=4,77 = 1.0 Gyr as a function of fimo and Hq. 
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FIG. 4. Contour map of zp which satisfies Atz=3M2 = 1.3 Gyr as a function of ^mo^^ and 
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FIG. 5. Contour map of zp which satisfies /S.tz=z.&2 = 1-5 Gyr as a function of ^moh^ and 
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FIG. 6. Contour map of zp which satisfies Ai2=4.77 = 1.0 Gyr as a function of Omo^^ and 
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